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The effect of a magnetic fieid on a gapped quantum magnet is described within the framework 
of the Quantum Dimer Model. A minimal model describing the proliferation of itinerant spinons 
above a critical field is proposed and investigated by Lanczos exact diagonalizations and quantum 
Monte Carlo simulations. For both square and triangular lattices, it is shown that spinons are fully 
polarized and Bose-condense. This offers a novel scenario of a Quantum Critical Point in the dimer- 
liquid phase (triangular lattice) characterized by the continuous appearance of a spinon superfluid 
density, contrasting with the usual triplet condensation picture. The possible role of other spinon 
kinetic terms neglected in the model are discussed. 
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Quantum (spin) disordered magnets are currently a 
very active field of research in Condensed Matter Physics. 
They are characterized by the absence of magnetic or- 
dering down to very low temperatures and, for some of 
them of interest here, by a spin gap to triplet excita- 
tions. In this context, the magnetic field is a particularly 
relevant external parameter as it can drive the system 
to a Quantum Critical Point (QCP) at which the spin 
gap vanishes [l[. In the conventional picture, the QCP 
can be viewed as the Bose condensation of triplets. Nu- 
merous examples of such QCP can be found in quasi-one 
dimensional magnets like the CuHpCl spin ladder Q, 
in two dimensional layered frustrated magnets like e.g. 
the Copper-Borate system Q or in the interesting three 
dimensional BaCuSiO dimer compound 

The resonating valence bond (RVB) liquid describes 
the paradigm for magnetically disordered systems. In- 
deed, RVB states describe pure quantum phases, which 
have no classical counterparts and, therefore, cannot be 
adiabatically connected to any band theory, based on 
Hartree-Fock or more sophisticated density functional 
approaches. Moreover, as suggested by Anderson in a 
cornerstone paper [Bj|, an RVB insulator could be inti- 
mately connected to high-temperature superconductors. 
In this field, an important step forward was done by 
Rokhsar and Kivelson, who formulated the RVB con- 
cept by introducing a simple microscopic Hamiltonian, 
the so-called Quantum Dimer Model (QDM) [|| . Simpli- 
fying the original electronic problem, they proposed an 
effective theory describing the dynamics of short-range 
dimers, that become the fundamental entities of the prob- 
lem. The QDM relies on two approximations which are 
believed not to alter significantly its physical relevance 
(and could be relieved at the price of bringing an enor- 
mous extra complexity): (i) only nearest-neighbor sin- 
glets are considered and (ii) different dimer configura- 
tions are taken to be orthogonal to each other, hence re- 



ducing singlets to "classical" dimers. However, quantum 
mechanics still enters in an essential way in the possibil- 
ity to flip plaquettes with parallel dimers. On the square 
lattice, the ground state of the QDM always shows lat- 
tice symmetry breaking. This fact has encouraged the 
(wrong) surmise that "fully" disordered ground state (in 
the sense of a state with neither SU(2) nor lattice symme- 
try breaking) cannot be stabilized in this framework. In 
this respect, the discovery that a gapped and disordered 
state can be obtained on the triangular lattice opened the 
way to an enormous number of investigations 0, Q . Re- 
cently, the zero-temperature properties of the QDM on 
the triangular and square lattices have been considered 
in great detail by numerical calculations [^, [l(| ■ 

In the QDM, both charge and spin excitations are 
frozen and the elementary excitations have topological 
lil 13j |. Some effort has been done in or- 
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der to introduce charged spinlcss excitations in the doped 
QDM, the so-called holons. The exact nature of these ob- 

H, O- Much more 
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jects has been widely debated 
subtle is the problem of inserting spinons, namely charge- 
less spin-1/2 objects, since, in this case, one faces a se- 

In fact SU(2) 



11] 



rious problem of non-orthogonality 
valence-bond configurations with unpaired spins at dif- 
ferent locations are not orthogonal to one another, in 
contrast to the holon case. The magnetic field offers a 
formidable tool to experimentally investigate such issues 
since it plays the role of a chemical potential w.r.t. the 
spinon density. 

In this Letter, we address the issue of spinon doping in 
both (lattice) ordered and disordered dimer backgrounds. 
Physically, this is achieved by the action of a magnetic 
field. On a general basis, we first derive all types of terms 
that should govern the spinon dynamics. Then, we pro- 
pose a simplified model where the (up or down) spinon 
can hop along plaquette diagonals (while its neighboring 
dimer rotates by 90 degrees) and where the dimer dynam- 
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FIG. 1: (Color online) The four different kinetic processes 
of the effective Hamiltonian and their corresponding matrix 
elements in the case of a (non-frustrated) Heisenberg quantum 
magnet, a = l/v2 is the singlet normalization factor and J 
the spin coupling. The sign convention for dimers is the one 
of Ref. H 
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FIG. 2: (Color online) The lowest-energy levels for four and 
six spinons on a 4 x 4 square lattice in the subspace of zero 
momentum. The couplings are t/J = 0.2 and v/J — 0.5 and 
states are classified according to the total spin S. 



ics is that of the standard QDM. By using Lanczos exact 
diagonalizations, we show that, within this model, for a 
finite magnetic field, all spinons are fully polarized, the 
energy gain being of kinetic origin as in the well known 
Nagaoka effect [Iff. The magnetization profiles are then 
computed using Green's function Monte Carlo. We find 
a new type of QCP characterized by Bose-Condensation 
of spinons in contrast to the usual scenario of triplet con- 
densate. The behavior of the magnetization at the QCP 
depends on the nature of the (spin) disordered ground 
state at zero field. 

Creating a pair of spinons in a gapped quantum mag- 
net (described by a QDM) will cost a finite energy (called 
Ab later on). However, a finite density x of (partially po- 
larized) spinons could be stabilized by the Zccman energy 
-gfXBH^Sf = — §iV(n T - n ; ), n T (n;) being the den- 
sities of up (down) spinons in the system, x = rij + 
and N is the number of sites. Therefore, as a preliminary 
work, we need to establish the form of the kinetic pro- 
cesses for the spinons moving in the fluctuating dimer 
background. At the microscopic level, such terms will 
a priori arise from the non-orthogonality of the basis 
set. Since the derivation for a true frustrated microscopic 
model is beyond the scope of this work, we simply con- 
sider here the (non-frustrated) Heisenberg model on the 
square lattice to get insights on the generic form of these 
matrix elements 17] . Notice that a similar type of deriva- 



tion has been curried out for a realistic but more complex 
model in a different context [l8j |. 

Following the same procedure as for the derivation 
of the QDM in the absence of spinons fl4j |. we calcu- 
late the largest contributions of the matrix elements, 
H c c i = (c' | TJncis | c) , where |c) and \d) are typical dimer 
configurations, including spinons. A careful analysis 
gives four elementary kinetic processes and their over- 
lap weights, as displayed in Fig. [TJ Interestingly, the 
processes involving a spin-flip of two spinons [term 3)] or 
two-spinon creation or annihilation [term 4)] have both 
positive and negative signs depending on the orientation 



of the spinons. This introduces a serious sign problem 
into the quantum Monte Carlo approach. In addition, 
the term 4), in contrast to the holon-doped case, does 
not conserve the number of "dopants" . 

Hereafter, we consider a simplified version of the QDM- 
spinon model, by retaining only the dimer flips (as in the 
usual QDM) and the spinon hopping 2). This model has 
only negative off-diagonal matrix elements and should 
be tractable by appropriate quantum Monte Carlo. In 
a second step, we shall discuss the possible role of the 
terms we have left behind. Note that in the case of a 
frustrated magnet, we expect the same type of processes 
although their magnitudes should depend on the precise 
microscopic magnetic interactions. Therefore, the QDM 
Hamiltonian is defined here by 

W =«£iV c |c><c|- l c ')(c|-i E l c ">H (!) 
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The sum on (c) runs over all configurations of the Hilbcrt 
space, N c is the number of flippable plaquettes, the sum 
on (c', c) runs over all configurations \c) and \d) that dif- 
fer by a single plaquette dimer flip, and the sum on (c", c) 
runs over all configurations \c) and |c") that differ by a 
single spinon hopping between nearest neighbors (trian- 
gular lattice) or next-nearest neighbors (square lattice), 
see Fig. [TJ The coupling to the external magnetic field is 
also added to the above Hamiltonian. 

To get insights on the ground-state properties of this 
system, we first exactly diagonalized (by standard rou- 
tines or a Lanczos algorithm) a 4 x 4 square cluster with 
periodic boundary conditions and with 4 and 6 doped 
spinons. Eigen-energies are labeled according to the to- 
tal spin quantum number (lij . Results are reported in 
Fig. [2] for zero magnetic field (h = 0). The important 
point is that, in all the cases we considered, the ground 
state has the maximum spin value, namely it is fully po- 
larized. This bears stron g si milarities with conventional 
Nagaoka fcrromagnetism [lfj (but for a finite concentra- 
tion of dopants, which is remarkable) in the sense that the 
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energy gain w.r.t. the singlet subspacc is of kinetic origin. 
Obviously, the fully polarized state is further stabilized 
by a finite magnetic field so that this property should 
be valid for all fields and dopant concentration making 
quantum Monte Carlo simulations available. Note how- 
ever that the neglected processes 3) and 4) of Fig.[T]might 
change the relative stability of the fully polarized state 
w.r.t. some lower spin states and this will be discussed 
later on. In this way, technically, the spin variable does 
not play any role and, therefore, one can borrow here the 
results obtained in the context of the monomer-doped 
QDM 2(], 21 ] , but accounting for the extra term involv- 
ing the magnetic field —hNx/2 (n-j- = x). It is of inter- 
est to notice that, in contrast to the case of hole-doping 
where Fermi statistics plays a role [Hj], "spinon-doping" 
do not introduce any minus-sign problem (here J > 0). 
In other words, since an electron destruction operator can 
be written as Cj CT = fjbi a in term of a fermionic holon- 
creation operator fj and a bosonic spinon-destruction op- 
erator bio-, the physical holon-doped and "spinon-doped" 
QDM simply differ from the statistic of the bare doped 
monomers. 

Here, we aim to assess the number of spinons that are 
present in the ground state, once an external magnetic 
field h is applied and competes with the dimer binding 
energy Ab- Indeed, —h plays the role of a chemical po- 
tential and controls the spinon density which gives the 
reduced magnetization. It is natural to introduce the ef- 
fective energy cost A°g = As — h of breaking a dimer 
into two unbound static spinons under a finite magnetic 
field. Note that the spontaneous appearance of spinons 
in the ground-state (which defines the QCP) should ap- 
pear for fields lower than the characteristic field h = Ab 
since spinons acquire extra kinetic energy compared to 
dimers (in other words, the actual critical value of A^f 
is positive). This study can be done by considering the 
energy per site e(x) of the holon-doped model [2l|, from 
which the Legendre transform F(x) = e(x) + ^A Q ^x can 
be constructed. The actual e(x) has been obtained by 
using the Green's function Monte Carlo method [22j], for 
both the square and the triangular lattices 2l| . The ac- 
tual spinon concentration xq for a given value of A^ is 
obtained by minimizing F(x) with respect to x. Notice 
that xq is nothing else but the ratio of the magnetizations 
M /A/sat, with M sa t the value at saturation, i.e., when the 
whole system is filled by spinons. 

The fact that the monomer-doped model may or may 
not show phase separation (PS), opens the way to two 
different scenarios for the behavior of the system as 
a function of the magnetic field. Let us first recall 
that, on the square lattice, as soon as v/J < 1, the 
ground state is unstable in the presence of an infinites- 
imal monomer concentration and phase-separates be- 
tween undoped (macroscopic) domains and regions with 
a concentration x c of monomers. Instead, in the trian- 
gular lattice, the infinitcsimally doped system is stable 
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FIG. 3: (Color online) Phase diagram of the QDM given by 
Eq. JT} as a function of magnetic field, for the square lat- 
tice (full squares) and the triangular lattice (empty triangles) . 
The dashed and continuous lines indicate first- and second- 
order transitions, respectively. 



for v/J> 0.9, whereas it phase-separates for smaller val- 
ues of the dimer-dimer repulsion [2l|. Such different be- 
haviors are reflected in the phase diagram of Fig. [3] in 
the (h — As,i) plane by continuous or first-order transi- 
tions between the regions with and without magnetiza- 
tion. Remarkably, the transition curves t = t c (h) have 
linear behaviors and, for a given lattice structure, the 
slope docs not depend upon the ratio v/J. Three differ- 
ent values of the dimer-dimer repulsion v/J have been 
chosen for the triangular lattice, two of them inside the 
stability region of the RVB liquid of the undoped system 
(v/J = 1 and 0.95), whereas the latter one (v/J = 0.85) 
corresponds to a region where the doped ground state 
phase separates, due to local dimer ordering. 

The behavior of the magnetization vs magnetic field 
can be understood on simple grounds. In a stable sys- 
tem, the compressibility is finite, implying that the sec- 
ond derivative of e(x) is positive (i.e., e(x) is a concave 
up function of x). In this case, e(x) ~ e(0) + ax + bx 2 /2, 
where b > from the stability criterion. Moreover, the 
linear coefficient is generally negative, i.e., a = —\a\. On 
the other hand, PS is signaled by an infinite compress- 
ibility and a vanishing second derivative of e(x). This 
leads to a precise form of the energy in the thermody- 
namic limit, namely e(x) = e(0) + ax for x < x c and 
e(x) ~ e(x c ) + a(x — x c ) + b(x — x c ) 2 /2 for x > x c . No- 
tice that the correct behavior of the stable case may be 
found from the previous one by imposing x c = 0. By 
minimizing F(x) we then find a general expression for 
the reduced magnetization, 
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for h > 2\a\ — Ab and M = otherwise. Therefore, in 
a stable system (i.e., x c = 0), by increasing the mag- 
netic field, there is a continuous transition between the 
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FIG. 4: (Color online) a) Spinon concentration for different 
values of t as a function of the magnetic field h (both in units 
of J) for the triangular lattice at v = 0.95. b) The same as 
in a), but for the square lattice and at v = 0.80. Dashed 
lines show the region where a first order transition occurs. 
Error-bars are of the size of the symbols. 



spin order. However, we believe that transverse magnetic 
ordering would not change our qualitative picture. 

In summary, a simple extension of the QDM including 
spinon doping is introduced to account for the Zeeman 
effect of an applied magnetic field. This model involves 
two semi-phenomenological parameters, the hopping am- 
plitude t for a spinon to hop (and exchange with a dimer) 
and the bare energy cost to break up a dimer into two 
spinous. Interestingly, the phase diagram of this model 
can be obtained from large-scale Green's function Monte 
Carlo methods jUJ showing, as a function of the field, 
either a first order transition (square lattice) or a new 
type of Quantum Critical Point (triangular lattice) char- 
acterized by Bose-Condensation of spinons. 
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gapped phase to a phase with a finite spinon density. In 
Fig. 2k., we show the behavior of x = M/M sat as a func- 
tion of h — Ab for the triangular lattice and v/J = 0.95, 
namely starting from a h = RVB liquid phase. A 
Bose-condensation of spinons is found, characterized by 
a second-order phase transition for all the values of t/J 
considered, hence confirming the analytical arguments. 

On the contrary, whenever x c ^ 0, by increasing the 
magnetic field, the number of spinons suddenly jumps to 
a finite value for h > Ab — 2|a|. In this case, the critical 
value of the corresponding (effective) dimer binding en- 
ergy is given by A^f c = 2|a|, i.e., twice the linear slope of 
the energy per site at zero doping, |a| = [e{x c ) — e(0)]/x c . 
In panel Fig. [4}), we consider the case of the square lat- 
tice, with v/J= 0.8. At small enough magnetic field, the 
ground state is an anisotropic mixed columnar/plaqucttc 
phase, as recently pointed out Upon increasing h, 

the ground state exhibits a first-order transition at which 
M/M sa t jumps suddenly to a finite value. 

Though serious ergodicity problems prevent us to con- 
sider very small values of t/J, we have evidence that, ex- 
cept for v/J= 1, the transition lines in Fig.[3]extrapolate 
to non-zero values of t/J as h — > Ab. This fact implies 
that the linear coefficient of the energy per site e(x) can 
be written as \a\ ~ at — /3(v/J). At the Rokhsar-Kivelson 
point v/J=l and for t = 0, we have that e(x) = 0, which 
implies (3(1) = and the vanishing of t c for h — > Ab- 
Away from the Rokhsar-Kivelson point, (3(v/J) is no 
longer vanishing and, from the fact that Aff c = 2\a\ (see 
above), we obtain t c (h) ~ (Ab — h)/(2a) + (3/a, that 
explains the behavior of the curves in Fig. [3] 

Lastly, we note that the neglected terms 3) and 4) 
(as well as a direct exchange of two neighboring up and 
down spinons) may produce a finite density of minority 
down spinons which in turn might induce transverse spin- 
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